Isobaric Yield Ratio Difference in Heavy-ion Collisions, and Comparison to Isoscaling 
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An isobaric yield ratio difference (IBD) method is proposed to study the ratio of the difference 
between the chemical potential of neutron and proton to temperature (A/^/T) in heavy-ion collisions. 
The Afj,/T determined by the IBD method (IB-A^/T) is compared to the results of the isoscaling 
method (IS-Ap/T), which uses the isotopic or the isotonic yield ratio. Similar distributions of 
the IB- and IS-Afi/T are found in the measured 140A MeV '"'''**^Ca + ^Be and the ^s.e^jyj _^ 9gg 
reactions. The IB- and IS-A/i/T both have a distribution with a plateau in the small mass fragments 
plus an increasing part in the fragments of relatively larger mass. The IB- and IS-Afi/T plateaus 
show dependence on the n/p ratio of the projectile. It is suggested that the height of the plateau is 
decided by the difference between the neutron density (p„) and the proton density (pp) distributions 
of the projectiles, and the width shows the overlapping volume of the projectiles in which pn and 
pp change very little. The difference between the IB- and IS-A^/T is explained by the isoscaling 
parameters being constrained by the many isotopes and isotones, while the IBD method only uses 
the yields of two isobars. It is suggested that the IB-A/i/T is more reasonable than the IS-A/i/T, 
especially when the isotopic or isotonic ratio disobeys the isoscaling. As to the question whether 
the A^/T depends on the density or the temperature, the density dependence is preferred since the 
low density can result in low temperature in the peripheral reactions. 

PACS numbers: 25.70.Pq, 21.65. Cd, 25.70.Mn 



I. INTRODUCTION 

In the models based on the free energy to predict the - 
fragments in heavy-ion collisions (HICs) above the Fermi • 
energy, the yield is mainly determined by the free energy, • 
the chemical potential of proton and neutron, the tem- • 
perature, etc. [H-ll]. In the ratios between the fragment • 
yields, some of the information which the fragment car- - 
ries will cancel out, and the retained information is useful - 
to study the properties of the colliding sources 0, [H , the ' 
fragment itself [y, 0, and the temperature of the reac- 
tions [ll,[3|- The isoscaling method is one of the important 
methods to constrain the symmetry energy of the nuclear ' 
matter in HICs [l, |3j 1 which makes it important for the 
study of the nuclear equation of state • The isoscaling ' 
phenomena are systemically studied experimentally, and ' 
extensively examined in theories from dynamical models ' 
to statistical models [Tl| - [2]| , The effects of the sec- ' 
ondary decay, which significantly influence the results, ' 
are also investigated [ij, [22| - [25j . Besides the isoscaling 
method, the isobaric yield ratio is promoted to determine 
the symmetry energy of the fragments produced in HICs , 
in a modified Fisher model [gI. [26l - l29| . At the same time, , 
the isotopic ratio and the isobaric ratio are also used to , 
study the temperatures of the colliding sources [l|, l30l - [37j , 
or the heavy fragments in HICs [33, [s^ . Since both , 
the isoscaling method and the isobaric ratio method are , 
deduced in the framework of the free energy theories, and , 
they both relate the yield of fragments to the symmetry , 
energies of the colliding sources, it is important to com- , 
pare the nuclear symmetry energy determined by them. 
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In this article, the difference between the chemical po- 
tentials of the neutron and proton will be compared using 
the isoscaling method and the isobaric ratio method. In 
Sec. nil the isoscaling and the isobaric ratio difference 
methods will be deduced in the framework of the grand- 
canonical ensembles. In Sec. IIIIl the fragment yields in 
the UOA MeV ^^''^^Ca + ^Be and 58,64]s^j _^ ggg reactions 
will be analyzed using the isoscaling method and the iso- 
baric yield ratio (lYR) method, and the results will be 
compared. A summary will be presented in Sec. IIVI 



II. MODEL DESCRIPTION 

The isoscaling and the isobaric ratio difference (IBD) 
methods in the grand-canonical ensembles will be intro- 
duced briefly. In the grand-canonical limit, the yield of a 
fragment with mass A and neutron-excess I (I = N — Z) 
is given by [H, 113 

Y{A, I) = CA^exp{[F{A, I) + fi„N + ^ipZ]/T}, (1) 

where C is a constant. N and Z are the neutron and 
proton numbers, t is nonuniform in different reaction 
systems [4l|. /i„ and /ip are chemical potentials of the 
neutron and proton, respectively; F{A,I) is the free en- 
ergy of the cluster (fragment), and T is the temperature. 

The isoscaling method is as follows: for one fragment 
in two reactions of the same measurements, based on Eq. 
(dl), the yield ratio of the two reactions, i?2f (-/V, Z), can 
be defined as 0, |3] 

RiUN, Z) = Y2{N, Z)/Yi{N, Z) = C"exp(a7V + PZ) 

(2) 

where C is an overall normalization constant which orig- 
inates from the different reaction systems. /i„ and /ip 



2 



are assumed to change very slowly; a = A/i„/T with 
A/i„ = At«2-Mni, and^ = A/ip/T with Ajj.p = ^J,p2~^J.pl, 
which reflect the properties of the colliding sources. In 
the isotopic ratios, /3 cancels out and a can be fitted; and 
in the isotonic ratios, a cancels out and /3 can be fitted. 
a « — /? is found, and a can be related to the symmetry 
energy (C'sym) nuclear mass of the colliding source by 

I Cs^™ r/ Zi \2 _ I Z2 \2 
> A 



"[("f^)'^ — {^Y], or some similar relationships : 



9 iHlfili^.' 

10 The isobaric yield ratio (lYR) is used to study the 39 

11 symmetr y e nergy of the fragment at finite temperatures 40 

12 [a, i, [Ulii]. When using lYR to study the /i„ and Hp, « 

13 the analysis method should be reconstructed. Starting 42 

14 from Eq. ([1]), in one single reaction, the lYR between 43 

15 the isobars differing by 2 units in /, R^^{I + 2, 1, A), can 44 

16 be defined as 45 

i?^^(/+2,/,A) = r(A/+2)/y(A,/) 

= exp{[F(/ + 2, A) - F{I, A) + - fip]/T}, (3) ^ 

17 The CA^ term in Eq. ([T]) cancels out and the system 

IS dependence is removed. Assuming that the isobars in ^° 

19 the ratio have the same temperature, only the retained " 

20 /i„ and lip are related to the colliding sources. Taking " 

21 the logarithm of Eq. one can obtain, 



InR^^il + 2, /, A) = {AF + Afi)/T, 



(4) 



where AF = F{I + 2,A)-F{I,A), and A^i = Hn-fJ-p- In i 
two reactions of the same measurements, the difference . 
between the lYRs, i.e. the IBD method, can be defined , 
as, ( 



Alni?^f 



= Afi^/T- 



-2,/, A)] 
Afip/T - 



-ln[R{^{I + 2J,A)] 
A/i/T = a - /?, (5) 



Eq. ^ also shows the relationship between the results of 
the isoscaling parameters {a and /3) and the IBD method. 
For convenience, the IBD and isoscaling A/i/T are la- 
beled as IB-A/(/r(= Alni?|f ) and IS-A/VT(= a - /?), 
respectively. 
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FIG. 1. (Color online) (a) The isotopic yield ratio of the 
isotopes from Z = 6 to 19, and (b) the isotonic yield ratio of 
the isotones from iV = 8 to 21 in the UOA MeV ""'"'^Ca + ^ 
^Be reactions The lines are the linear fitting results. 



III. RESULTS AND DISCUSSIONS 

The yields of the fragments produced in the 140A MeV 
^O'^^Ca + ^Be and the ^^'^'^Ni + ^Be reactions were mea- 
sured by Mocko et al. at the National Superconducting 
Cyclotron Laboratory (NSCL) in Michigan State Uni- 
versity. The details of the measurements were described 
in Ref. [i^. The isoscaling phenomena in these reac- 
tions were studied in Ref. [IJ], and the isotopic (iso- 
tonic) yield distributions in these reactions were stud- 
ied in Ref. [i^. In this article, the IB- and the IS- 
A^i/T associated with the fragments in these reactions 
will be analyzed. The analysis will be performed be- 
tween the isotopic ■*®Ca/'*°Ca and ^^Ni/^^Ni reactions, 
the n/p symmetric ^*Ni/*°Ca reactions, and the neutron- 
rich ^^Ca/^'*Ni reactions. The reaction of the relatively 
small n/p projectile is denoted as 1, and the other one 
as 2. The isoscaling parameters a and /? are obtained 
from the linear fitting of the isotopic ratio and the iso- 
tonic ratio in the chosen reactions according to Eq. ([2|) . 
For example, in Fig. [1] the isoscaling phenomena of the 
fragments in the lAOA MeV ^°'^^Ca + '^Be reactions are 
shown. In Figs. [IJa) and (b), the isotopic scaling and 
the isotonic scaling are plotted, respectively (similar re- 
sults can be found in Ref. [13|). a {(3) equal to the 
slope of the linear fitting of the isotopic (isotonic) scal- 
ing. For one fragment, the IS- A/i/T is calculated using 
the a and /? obtained from its Z isotopes and N isotones, 
respectively. For IB-A/i/T, the lYR in each reaction is 
calculated first, then the difference between the lYRs in 
the two reactions is calculated according to Eq. 
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FIG. 2. (Color online) The IB- and IS-A/i/T in the 140A MeV 
*°''**Ca -f ^Be reactions I = N — Z is the neutron-excess. 
The lines are for guiding the eyes to the plateaus. 



In Fig. m the IB- and IS-A/i/T in the ^^Ca/^^Ca + 
^Be reactions are plotted. Very similar trends of the IB- 
and IS-A/t/T distributions in each /-chain of fragments 
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are found. In each /-chain, both the IB- and IS-A/i/T in 
the small-^ fragments form plateaus (around A^/T = 
2). and Afi/T increases as A increases in the large- A 
fragments. The plateuas of the IB- and IS-A/i/T almost 
overlap, and it is interesting that the IB- and IS-A/i/T 
have very little difference in the / = [panel (b)] and 
the 1=1 [panel (c)] fragments. Except for the / = 
and I ~ I fragments, differences between the IB- and 
IS-A/i/T in the large-^ fragments are shown. 
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FIG. 3. (Color online) The IB- 
MeV ^ 9gg reactions [42 

the eyes to the plateaus. 



and IS-A/i/T in the UOA ' 
. The lines are for guiding ' 



In Fig. m the IB- and IS-A/i/T in the ^^Ni/s^Ni + ' 
^Be reactions are plotted. Very similar results as the ' 
^®Ca/^°Ca reactions are found, except that the values of' 
the plateaus decrease to about 1.4. The values of the IB- ' 
and IS-A^/T almost overlap in the 7=1 [panel (b)] and ' 

15 the 1 = 2 [panel (c)] fragments. The n/p for *®Ca/^°Ca ' 

16 is 1.4/1.0, which is larger than the value 1.286/1.071 for < 
64Ni/58Ni. Comparing to the results in the '^^Ca./^^Ca- 
reactions, obviously larger widths of the IB- and IS-A/i/T ' 
plateaus in the ^*Ni/^®Ni reactions are found. ' 

From the results shown in Figs. [5] and [31 for the ' 
plateau, its height (the value) and its width (the nu- ' 
clei ran ge i t covers) should be noticed. In the statistical ! 
models |43l . [4^ , the yield of a fragment to some extent is ! 
decided by the density distributions of protons (pp) and ! 
neutrons (p„) of the projectile and the target nuclei. A: 
nucleus can be assumed to have a core region, in which ! 
Pp and pn change very little, and a skirt region, in which e 
Pp and Pn change fast. The Pp distributions of isotopes e 
can be assumed to be similar, especially when the masses e 
of the isotopes do not differ much. This indicates that, e 
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FIG. 4. (Color online) The IB- and IS-A/i/T in the UOA 
MeV '"'Ca + ^Be and '^^Ni + ^Be reactions The lines 

are just for guiding the eyes of the plateaus. 



the height of the plateau is decided by the difference be- 
tween Pn and Pp in the projectiles, and the width shows 
the overlapping volume of the projectiles in which p„ and 
Pp vary slowly. Regarding the A/i/T in the isotopic pro- 
jectile reactions, the height of the plateau indicates the 
difference between the pn distributions of the projectiles. 

In addition to the similarity of the isotopic distribu- 
tions in the '^°Ca/''*Ca and the ^^Ni/^^^Ni reactions, the 
similarity of the isotopic and the isotonic distributions 
in the ^^Ni/^°Ca and the ^®Ca/^*Ni reactions were also 
investigated in Refs. [i^, [4^. It is found that the iso- 
topic or isotonic yield distribution shows dependence on 
Pn and Pp of the projectiles. In the isoscaling method, the 
system effects do influence a and /3. In the IBD method, 
the system effects are also removed in the lYR. Thus the 
A/t/T in the n/p symmetric ^^Ni/'^'Ca reactions, and the 
neutron- rich '^^Ca/^^Ni reactions can also be analyzed. 
In Fig. [H the A/i/T in the ^*Ni/''°Ca reactions are plot- 
ted. The n/p of 58Ni/40Ca is 1.071/1.0, thus /)„ and pp 
for them can be assumed to have similar trends but differ 
in values. In Fig. [31 very small differences between the 
IB- and IS-A/i/T are found in each Tchain. The values 
of the plateaus decrease to about 0.7, and in the / =0 
and 1 chain, the IB-plateaus become even smaller. 

The IB- and IS-A/i/T in the ^^^^Ca/^^Ni reactions are 
plotted in Fig. M The n/p ratio of ^^Ca./^^m is 
1.4/1.286. Though large differences between the IB- and 
IS-A/i/T in the I < 2 fragments are shown, similar values 
of the IB- and IS-A/i/T in the / > 3 (A > 25) fragments 
are found. The plateaus in the / < 3 chains decrease 
to smaller than 0.5. The characteristics of the IB- and 
IS- A/t/T distributions are very similar to those of the 
isotopic or isotonic distributions shown in Rcf. [43[, i.e.. 



4 



the density dependence of p„ and pp in the projectiles. 
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FIG. 5. (Color online) The IB- and IS-A/i/T in the 140A 
MeV ''*Ca + ''Be and "^^Ni + ''Be reactions [H. The hnes ' 
are for guiding the eyes to the plateaus. 



In Fig. [2]ini similarities between IB- and IS-A/i/T are . 
shown, i.e., the plateaus and the values in the small- . 
A fragments of the ^^C&/^°Ca. and ^^^Ni/^SNi reactions, ' 
most fragments of the ^^Ni/^°Ca, and large-A fragments ^ 
of the ^^Ca/^^'Ni reactions; while the differences between ^ 
the IB- and IS-A^/T are shown in the large-^ fragments ^ 
of the 48Ca/40Ca and ^^^Ni/^^Ni reactions. The values' 
of the plateaus show a dependence of the n/p ratio of 
the reaction systems, which are about 2, 1.4, 0.7, and 0.5 ' 
in the 48Ca/40Ca, 64Ni/58Ni, ^^Ni/^Ca, and "SCa/^^Ni ' 
reactions. ' 
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FIG. 6. (Color online) The IB- and IS-A/i/T between the , 
fragments in the impact parameter regions of Rl (6=0-2 fm) , 
and R2 (b=6-8 fm) in the 140A MeV ^^Ki + ''Be reactions , 
calculated using the AMD -I- GEMINI models. 



To see more clearly on this point, we performed a sim- 
ulation of the 140A MeV ^^Ni -|- ^Be reaction using a 
microscopic transport model. Though there are many 
choices- such as the quantum molecular d yna mics (QMD 
and its improved versions) models [2l|, [22I.I33. |47| - |49| , the 
AMD model [H, El HH, [12] , and different methods to 
form the cluster or fragment [13, the AMD model 
plus the sequential decay GEMINI code [HI] have been 
used to simulate the reaction, since similar works are 
performed, and the experimental yields of fragments are 
well reproduced [i^. In the calculation, the standard 
Gogny interaction (Gogny-gO) is used [s^l, the fragments 
are formed using a coalescence radius, Rc = 5fm in the 
phase space at the time t — 500/m/c in AMD. Two cuts 
of impact parameters arc used in the fragments analysis, 
i.e., 6=0-2 fm (labeled as Rl) and 6=6-8/m (labeled as 
R2). The results are plotted in Fig. [6] Since ^^Ni is a 
symmetric nucleus, its p„ and pp distributions are very 
similar, and the densities only decrease sharply in the 
very edge of the nucleus. In Rl and R2, the difference 
between p„ and pp is also very small. According to the 
assumption that the plateau of A/i/T depends on n/p 
of the projectiles, the plateaus of the IB-Ap/T should 
be very small. In Fig. [6l it can be seen that the IB- 
and IS-Ap/T are very similar, and their values are very 
small. The IS-A/i/T of most fragments are in a range 
of 0.1 ± 0.1, and the IB- A/i/T of the / = - 2 chains 
arc 0± 0.1. This indicates that the IB- and the IS-A/i/T 
have very little difference when pn and pp of the two pro- 
jectiles are similar, and differ if /?„ and pp of the two 
projectiles are different. 

Furthermore, the difference between the IB- and IS- 
A/t/T should also be discussed. Generally, the IB- and 
IS-A/t/T should be the same since they are obtained from 
the same fragments and in the same theory. Though 
in most fragments the IB- and IS-A/z/T are very sim- 
ilar, the difference between the IB- and IS-A/t/T are 
also shown in many fragments. For examples, in frag- 
ments where the IB- and IS-A/i/T are different, the IB- 
A/t/T are larger than the IS-A/i/T in the 4*^Ca/4°Ca and 
^'^Ni/^^Ni reactions, while the IB-A/i/T arc smaller than 
the IS-A/i/T in the ^^Wi/^^Ca. and ^SCa/^^Ni reactions. 
According to Eq. ([2]), the isoscaling parameters a (/3) 
are obtained from the linear correlation between the iso- 
topic (isotonic) ratio and neutron (proton) numbers, and 
a {(3) is the same for all the isotopes (isotones) . In other 
words, a {(3) is the scaled parameter for all the isotopes 
(isotones). For one fragment, its a (/3) is constrained 
by its isotopes (isotones), thus a (/?) can not reflect the 
difference between the isotopes or isotones. While ac- 
cording to Eq. (I5|), the IB-A/i/T result only relies on the 
two related isobars, and the difference between isobars of 
different masses can be obviously shown. It is suggested 
that since the IBD method uses only two isobars, the 
IB-A/i/T result is not influenced by the rest fragments, 
and more precise results than the isoscaling method can 
be obtained, especially when the fragment disobeys the 
isoscaling. 



Finally, we discuss the temperature effect in the Afi/T. so 
Hn and fj,p depend on both the density and the temper- 51 
ature. In theories based on the free energy, it is difficult 52 
to separate the free energy and the temperature 0, [l^l ■ 
In the IB- and the IS-A/i/T, it is also difficult to sep- 
arate A/.i and T. Besides the density effects in A/i/T, 
the temperature also influences A/i/T. Actually, the 
temperature should be defined at thermal equilibrium, 
but in intermediate energy HICs no thermal equilibrium 54 
is reached. In other words, the temperature IS nonuni- 55 
form in the collisions. In QMD, the "temperature" cause 
be extended to the non-equilibrium situations and ex- 57 
tracted in the local density approximation [s^ . In grand- 58 
canonical ensemble theory, the temperature is supposed 59 
to be the same, but differs in each reaction system. In eo 
a recent work using a canonical thermodynamic model, ei 
a temperature profile of impact parameter (b) is intro- 62 
duced, in which the temperature decreases as b increases 63 
[5^. Considering the multiple sources collisions of dif- m 
ferent strengths according to the impact parameters, the 65 
temperature changes with the excitation energy. In the m 
Fcrmi-gas relationship the correlation between the ex- 67 
citation energy per nucleon {E*/A) and temperature is 68 
E*/A = rVa, or r = y^E*a/A, in which a = Ak and k 69 
is the inverse level density parameter. In Ref. [13], in the 70 
'**Ca/'^°Ca -I- ^Be reactions, a is found to decreases when 71 
E*/A increases (which corresponds to T = 1.2 — 2.14 72 
MeV), but a tends to be similar if £'*/yl is high. Noting 73 
that a Pi — /?, the temperature dependence of a can also 74 
explain the plateau plus the increasing part of the IB- 75 
and IS-A/i/r distributions as follows: if A/i are uniform 76 
in the source, the plateau forms in the fragments which 
have high E* /A, while the A^jl/T increases when E* /A 
decreases in the fragments which have low E* /A. In the 
statistical abrasion-ablation (SAA) model, the excitation " 
energy is E* = 13. 3A^ MeV, in which A^ is the number 
of nucleons removed from the projectile by the ablation- 78 
abrasion process 0, [H]. Then T = y/U.ikAA/ A. In 79 
peripheral collisions, due to the low density of nucleons, 80 
the abraded nucleons are fewer than those in the central 81 
collisions, which results in the relative low temperature. 82 
To some extent, the low temperature is a result of the 83 
low density in the peripheral collisions, which is similar 84 
to the temperature profile in Ref. [56| . Thus low density 85 
can also result in an increase of A/x/T. In an isobaric 86 
method, the temperature of the measured heavy frag- 87 
ment is suggested to be similar due to the significant in- 88 
fiuence of the secondary decay process Q. Thus, though 89 
the density dependence and the temperature dependence 90 



of A/x/T cannot be totally separated, the density depen- 
dence is preferred since the low temperature is one result 
of the low density. 



IV. SUMMARY 

To summarize, in the article, a new IBD method is 
proposed to investigate A/i/T of the colliding sources, 
and the result is compared with the result of the usually 
used isoscaling method within the same grand-canonical 
model. The IB- and IS-A/z/T are found to be similar 
in the distributions, which both have a plateaus in small 
mass fragments plus an increasing part in relatively larger 
mass fragments. The IB- and IS-A/i/T plateaus show 
dependence on the n/p ratio of the projectiles. It is sug- 
gested that the height of the plateau is decided by the 
difference between /)„ and pp in the projectiles, and the 
width shows the overlapping volume of the projectiles in 
which pn and pp change very little. The difference be- 
tween the IB- and IS-A/i/T is explained by a and /? being 
constrained by the many isotopes and isotones, while the 
IBD method only uses the yields of two isobars. It is 
suggested that the IB- A/i/T is more reasonable than the 
IS-A/i/T, especially when the isotopic or isotonic ratio 
disobeys the scaling. As to the question whether A/i/T 
depends on the density or the temperature of the col- 
liding, the density dependence is preferred since the low 
density can result in the low temperature in the periph- 
eral reactions. 
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